Controlling edge states in the Kane-Mele model via edge chirality 
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We investigate the dependence of band dispersion of the quantum spin Hall effect (QSHE) edge states in 
the Kane-Mele model on crystallographic orientation of the edges. Band structures of the one-dimensional 
honeycomb lattice ribbons show the presence of the QSHE edge states at all orientations of the edges given 
sufficiently strong spin-orbit interactions. We find that the Fermi velocities of the QSHE edge-state bands 
increase monotonically when the edge orientation changes from zigzag (chirality angle 6 = 0°) to armchair 
(9 — 30°). We propose a simple analytical model to explain the numerical results. 



In their seminal paper fit], Kane and Mele proposed a sim- 
ple two-dimensional model which realizes the quantum spin 
Hall effect (QSHE). The model essentially considers a tight- 
binding model on a honeycomb lattice akin to graphene with 
added spin-orbit interactions. Realistic graphene is charac- 
terized by only very weak intrinsic spin-orbit coupling of the 
order of 10~ 5 eV ||2|-|6[]. While no QSHE has been experimen- 
tally observed in graphene, the Kane-Mele construction has 
become a popular model of Z 2 topological insulators 13, Si- 
Significant attention is currently devoted towards understand- 
ing the relationship between the crystallographic orientation 
of edges and surfaces in topological insulators and the re- 
sulting properties of topologically non-trivial boundary states 
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In this work, we establish a dependence between the band 
dispersion of the QSHE edge states in the Kane-Mele model 
and the crystallographic orientation of the edges. In particular, 
we show that Fermi velocities of the topological edge states 
strongly depend on edge orientation. 

We investigate the electronic band structures of ID peri- 
odic honeycomb lattice ribbons within the Kane-Mele model 
Hamiltonian (jfl 
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where and indicate first and second nearest 

neighbors, respectively, and a is the spin index. In this ex- 
pression, the first term corresponds to an ordinary nearest- 
neighbor tight-binding model with hopping energy t. The 
second term introduces spin-orbit coupling of strength t so , 
Vij = ±1 is the Haldane factor 111 ill defined as = (dik x 
djk)/\dik x djk\ for a pair of second nearest neighbor sites 
{{i,j)) connected via a common neighbor k. a z is a Pauli 
matrix describing electron spin. In graphene-like systems, the 
spin-orbit term opens a gap A so = 6\/3£so at the Dirac points. 
Following Kane and Mele [1], we choose the spin-orbit sec- 
ond neighbor hopping t so = 0.03i. We stress that the value 
used overestimates the intrinsic spin-orbit coupling present in 
realistic graphene, which was predicted to be between 1 and 
50 /ieV according to first-principles calculations JJ-lfll. 

The configurations of investigated ID periodic honeycomb 
lattice ribbons are defined by two parameters: (i) the crys- 
tallographic orientation of the edges and (ii) the width of the 
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FIG. 1: (a) Atomic structure model of a chiral honeycomb ribbon 
characterized by the (4,1) edge translation vector (9 = 10.9°) and 
width w = 3. (b) Schematic illustration showing the projection of 
points K and K' of the 2D Brillouin zone of honeycomb lattice onto 
the direction corresponding to a chiral edge. Schematic band struc- 
tures of (c) zigzag and chiral ribbons (0° < 9 < 30°), and (d) of 
armchair ribbons (9 = 30°) in the presence of spin-orbit interac- 
tions. The quantum spin Hall effect edge states are shown as solid 
lines. 



ribbon. The edge direction is described by a translation vector 
(n, m) of the graphene lattice (see Fig. QJ). The high sym- 
metry directions, armchair and zigzag, correspond to vectors 
(1, 1) and (1, 0), respectively. Equivalently, the edge orienta- 
tion can be described in terms of chirality angle defined as 
the angle between the edge and the zigzag direction 1 12 , . 
The edge translation vectors and chirality angles are related to 
each other by the following relation: 
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Following Ref. [ 12], we defined the width w of the ribbon by 
the vector (—w, 2w) along the armchair direction as shown on 
FigOJ. 

Figure 12 shows the band structures of honeycomb ribbons 
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FIG. 2: Sequence of band structures of the honeycomb lattice ribbons characterized by different periodicity vectors with widths defined by 
«; = 16 and t so = 0.03t. Band structures calculated in the absence of spin-orbit interactions are shown for the cases of zigzag (8 = 0°) and 
armchair (8 — 30°) nanoribbons. The scales of the plots account for the varying Brillouin-zone dimensions. 



with chirality angles ranging from 9 = 0° ((1, 0) zigzag edge) 
to 9 = 30° ((1,1) armchair edge) via a series of intermediate 
edge orientations (chiral edges). All considered models have 
comparable width defined by w = 16. In Figure |2j one can 
immediately notice that all band structures feature linear band 
crossings occurring either at k — or at the Brillouin-zone 
boundary k = n/a. The crossings display a clear increase 
of the Fermi velocity v-p upon increasing 9. This relationship 
will be discussed in detail below. Analysis of the electronic 
states at the band crossings reveals that the channels of op- 
posite spins are localized at the opposite edges of ID ribbon 
structures. That is, all investigated ID honeycomb ribbons 
are in quantum spin Hall phase and exhibit spin-filtered edge 
states topologically protected against backscattering. 

The effects of spin-orbit term are clearly illustrated for the 
case of a zigzag (1,0) edge shown in Fig. [2] In the ab- 
sence of spin-orbit coupling (t so = 0; red dashed line) the 
band structure of the ribbon model exhibits a dispersionless 
band at E = 0. This band is four times degenerate (2 spins 
x 2 edges); it corresponds to edge-localized states originat- 
ing from the lifted compensation between the two sublat- 
tices of the honeycomb lattice [14]. The flat band connects 
k = 27r/3ao and k = — 27r/3ao (ao is the lattice constant 
of the honeycomb lattice). These momenta correspond to the 
projections of points K and K' of the hexagonal Brillouin 
zone (the locations of the Dirac cones in the band structure of 
graphene) onto the momentum space of the ID ribbon struc- 
tures (points K and K' in Fig.Ht). The introduction of spin- 
orbit term opens a band gap A s0 at K and K' lifting the de- 
generacy of edge states and leading to a non-zero value of vf- 
In the quantum spin Hall phase the edge states connect the va- 
lence band at K with the conduction band K', and vice versa 

(Fig.nt). 

The increase of edge chirality angle 9 has a distinct effect 
on the electronic structure of honeycomb ribbons as it reduces 
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FIG. 3: Fermi velocity vp of the quantum spin Hall effect edge states 
as a function of chirality angle 8. The labels indicate the correspond- 
ing periodicity vectors of the edges. 



the separation between points K and K' (see Figs. [U),c). 
More precisely, the distance between points K and K' is given 
by EEli 
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A = sin(7r/6- 9). 
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This allows us to provide an estimate of the Fermi velocity as 
a function of spin-orbit interaction strength t so and chirality 
angle 9: 



vf 



(4) 
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Figure [3] compares the magnitudes of vf obtained from 
band structure calculations performed on w — 20 honeycomb 
ribbons, with the estimates provided by analytic expression 
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©. For 0° < 8 < 20°, the analytic formula © shows very 
good agreement with the numerical results. As the chirality 
angle approaches the armchair edge limit, the computed val- 
ues of vp deviate from analytic estimates eventually resulting 
in a finite Fermi velocity at 8 = 30°. The case of armchair 
edges is special as both K and K' are projected onto k = 
in the ID ribbon band structure. The edge state dispersion in 
this situation is illustrated in FigQJl as well as in the calcu- 
lated band structure in Fig. [2] It follows that the Fermi veloc- 
ity of the linear edge-state bands recovers the Fermi velocity 
of the massless Dirac fermions in the bulk when spin-orbit in- 
teractions are absent (the case of graphene), v F = 
(Fig. [3}, confirming the recent result of Gosalbez-Martmez 
et al. lfl7n . Interestingly, this result does not depend on 
the strength of spin-orbit interactions, contrary to the low- 
8 regime. On the other hand, armchair ribbons as well as 
high-0 chiral ribbons of finite width are semiconducting in the 
absence of spin-orbit interactions 112, 3. Thus, spin-orbit 
coupling above certain critical strength is required in order to 
bring these systems into the quantum spin Hall regime. 

In summary, we investigated the dependence of the band 
dispersion of the topologically non-trivial edge states in the 
Kane-Mele model on the crystallographic orientation of the 
edges. It was shown that the Fermi velocity of the quantum 
spin Hall edge states increases monotonically upon varying 
the edge chirality angle from 6 = 0° (zigzag edge) to 8 = 30° 
(armchair edge). A simple analytical model estimates the min- 
imum Fermi velocity as v-p = 9\/3t so ao/n. The maximum 
value achieved for armchair edges recovers the Fermi velocity 
Vp = \/3too/2 of the Dirac fermions on honeycomb lattice 
in the absence of spin-orbit interactions. The relations estab- 
lished for this prototypical topological insulator provide an 
important insight into tailoring the properties of topologically 
protected boundary states in realistic materials. 
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